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We use distribution amplitudes of the light T, baryon and the most general form 
of the interpolating current for heavy T,b baryon to investigate the semileptonic Sfe — )• 
S/"*"/" transition in light cone QCD sum rules. We calculate all twelve form factors 
responsible for this transition and use them to evaluate the branching ratio of the 
considered channel. The order of branching fraction shows that this channel can be 
detected at LHC. 
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1 Introduction 



The systems involving heavy quarks decays are important frameworks to restrict the stan- 
dard model (SM) parameters as well as search for new physics effects. Especially, the flavor 
changing neutral current (FCNC) transition of 6 — )■ which is underlying transition 
of Sf, — )■ E/"*"/" decay at the quark level, is known to be sensitive to new physics effects. 
This process can also be used in exact determination of the Vtb and Vts as elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and answering some fundamental questions 
such as CP violation. 

In the last decade, important experimental progress has been made in identification 
and spectroscopy of the heavy baryons with single heavy quark [1]-|8]. It is expected that 
the LHC will open new horizons not only in the identification and spectroscopy of these 
baryons, but also it will provide possibility to study the weak, strong and electromagnetic 
decays of heavy baryons. 

In accordance with this experimental progress, there is an increasing interest on calcula- 
tion of parameters of the heavy baryons and investigation of their decay modes theoretically. 
The masses of these baryons have been calculated using various methods such as quark mod- 
els P4T7]. heavy quark effective theory [18142^ and QCD sum rules [251433] . Besides the 
mass spectrum, their weak, strong and electromagnetic decays have also received special 
attention, recently (for instance see and references therein). 

In the present work, we analyze the semileptonic Sf, — )■ S/"*"/" transition in the framework 
of the light cone QCD sum rules. The main ingredients in analysis of this channel are 
form factors entering the transition matrix elements. Using the most general form of the 
interpolating field for the heavy baryon as well as the distribution amplitudes (DA's) of 
the light S baryon, we first calculate all twelve form factors in full theory. Then, we use 
these form factors to calculate the total decay rate as well as the branching ratio of the 
considered decay channel. 

The paper is organized in three sections. In the next section, we obtain QCD sum rules 
for the form factors. In section 3, we numerically analyze the form factors and use them to 
calculate the related decay rate and branching fraction. 

2 light cone QCD sum rules for form factors 

In this section, we focus on the calculation of the form factors corresponding to — )■ Hl^l^ 
semileptonic decay which proceeds via 6 — )■ s transition at quark level. The effective 
Hamiltonian describing this transition is written as: 
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The amplitude of the transition can be obtained by sandwiching the effective Hamiltonian 
between the initial and final states, 

M = - l^)m^)hn + C,, (S|.-7^(l - 75)6|S,)/7^5/ 

- 2m5 Ct^ (S|s^(7^,g^(l + ^^)h\L^)h^l^ . (2) 

From this equation, it is obvious that the transition matrix elements (S(p) |s7^(l— 75)6|Eb(p+ 
q) and (S(p)|scr^,^g'^(l + ^^h\L\,{j) + g)) are required. These matrix elements are expressed 
in terms of twelve form factors /j, gfj, jj and gj [i running from 1 to 3) as follows: 

(S(p) I ^7^,(1 - 75)& I + q)) = u^ip) [fMq') + ta^.q" f2iq') + q^fsiq') 

- lulbQiiq^) - icr,,ul5q''92{q'^) - q^l593{q'^)\ uj:^,ip + q) , (3) 

and 

(E(p) I sia,,q''{l + 75)& | + q)) = uj:{p) [f^f^q') + ia^.q^f^q') + q'fi(q') 

+ ini^gliq^) + ifj^^uim^'gliq^) + q^i^gliq^) «s,(p + q) , (4) 

where ms^, and u-£ are the spinors of T,b and S baryons, respectively, and q denotes trans- 
ferred momentum. 

Our main task is to calculate the aforesaid transition form factors. In accordance with 
the philosophy of QCD sum rules, we start considering the following correlation functions: 

ni(p,g) = ^ J d'xe-''^^{Q I T{J^^(0),6(x)7^(l - 75)^(^))} I , 

K'iP^<i) = ' J d'xe-"^^{Q I T{J^H0),6(a;)^a^.g^(l + 75)5(0;)} | E(p)) , (5) 

where J^*" stands for interpolating current of E{,. The interpolating current should be chosen 
as a composite operator that has the same quantum numbers as the baryon under study. 
For Ef, baryon, there are two possible choices for such a current that docs not contain any 
derivatives or auxiliary four vectors. The most general form for the interpolating current is 
a superposition of these two choices. Hence, for the interpolating current of the E;, baryon, 
the operator 

J^^(x) = ^e,,,|[<(x)C6^(x)]75rf'^(x)+^[<(x)C75&'(x)]rf'^(x) 

- [h'''^{x)Cd\x)]^^u%x) - l3[U'^{x)C^^d\x)]u^{x)^ , (6) 

is chosen. Here C is the charge conjugation operator, /3 is an arbitrary parameter, and a, 
6, and c, are the color indices. Taking /3 = — 1 corresponds to the loffe current. 



2 



The correlation function can be calculated both in terms of the hadronic parameters, 
such as the form factors, and also in terms of the QCD parameters. The expression in terms 
of the QCD parameters is evaluated by expanding the time ordered product of the currents 
in terms of the E distribution amplitudes via operator product expansion (OPE) in deep 
Euclidean region. On the other hand, the physical counterpart is calculated by inserting a 
complete set of intermediate states. The two expression are then matched using dispersion 
relations. 

To begin with, let us evaluate the correlation function in terms of hadronic parameters. 
After inserting the complete set of intermediate states into the correlation functions and 
isolating the ground state contribution we obtain 



[p, q> = y 5 7 — \ — — \ — , (8) 

where the ■ ■ ■ stands for contributions of the higher states and continuum, and the sum is 
over the polarizations of the Sj, baryon. The matrix element of the interpolating current 
between the vacuum and the E;, baryon appearing in Eqs. ([7]) and (E]), (0 | J^''(0) | + 
q,s)), can be expressed in terms of the residue of the T,b baryon defined as: 

(0 I J^''(0) I + g, s)) = As,ns,(p + g, . (9) 

The other matrix elements in Eqs. ([?]) and (|8]) are defined in terms of the form factors as 
previously shown. Combining Eqs. (E]), (jlj), and (l7])-(j9]) and summing over the polarization 
of the Sb baryon using the expression 

X^ME,(p + g,s)MSi,b + g,s) ^+ms, , (10) 

s 

the correlation functions can be expressed as: 

n^(p, q) = ^^^ ^2 _(p^g)2 pA^/i - h + qjs 

-it^i^di - icrf,uq''i592 + qf^i^gs^uM , (n) 

+ 7^75^^ + i(J^ul5q''92 - q^ll59^}uJ:{p) . (12) 

Commuting ]^ all the way to the right and using the equation of motion to write ^uy,{p) = 
msu-z^p), Eqs. (fTTj) and (fT2|) lead to the final expressions for the phenomenological side: 



ni(p,g) = ^2 y^p^^y {Viiq')P, + mq')P, 4+ [f2{q') + h{q 
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+ other structures \uy.{p) , 



f 75 



(13) 



"^1, - (P + 



9l{q')+9l{q') 



+ other structures I Me (p) 



(14) 



In these two expressions only the independent structures, p^, p^^, g^^, P/^^75, and 

g^^75, are presented exphcitly, owing to their sufficiency to determine the aimed form 
factors, /i(/f ), /2(/J), /2 + /3(/2^ + fl). giigl), 92i9l) and ^72 + 93(9! + 9l)- 

After completing the evaluation of the correlation function in terms of the hadronic 
parameters, now let us focus our attention on evaluating the correlation function in terms 
of the QCD parameters and the DA's of the E baryon . After placing the explicit expression 
of interpolating current given in Eq. ([6]) into Eq. ([5]) and contracting out the heavy quark 
operators, we attain the following representation of the correlators in QCD side: 



n 



—I 



V2 



(C)r,/3(75)p0 - (C)w(75; 



+ /3 
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II 
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J cry ' 



The heavy quark propagator, Sq{x) is calculated in 



(16) 



where. 



[X) - IQs 
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and Ki are the Bessel functions. Note that Sq'^'^ represents the free propagation of the 
heavy quark, and the remaining terms represent the interaction of the heavy quark with 
the external gluon field. The calculation of the contributions of the latter effects require 
the four- and five-particle baryons DA's which are currently unknown. But since they are 
higher order contributions, they are expected to be small [15H17] and we ignore them in 
the present work. In [H], it is also found that the form factors entering the semileptonic 
decays of the heavy A baryons turn out to receive only a very small contribution from the 
gluon condensate. 

The matrix element e"^^(0|M^(0)s^(a;)(i^(0)|E(p)) can be expressed in terms of S baryon's 
wave functions and are given in [IH], and for completeness explicit form of them are pre- 
sented in the Appendix. After evaluating the Fourier transform, the correlation function is 
expressed in terms of the QCD parameters and the DA's of the E. 

The sum rules are obtained by first Borel transforming both expression of the correlation 
functions and then equating the coefficient of various structures. Finally, the contributions 
of the higher states and the continuum are subtracted using quark hadron duality. To 
extract the numerical value of the form factors, value of the residue is also required. The 
residue of the J2f, baryon is calculated in [50] . 

3 Numerical Analysis 

In this section, we perform numerical analysis of the form factors and use them to pre- 
dict the decay rate and the branching ratio. The masses of the Sf,, S baryons and the b 
quark are taken as ms, = (5807.8 ± 2.7) MeV |52], ms = (1192.642 ± 0.024) MeV, and 
mb = (4.7 ± 0.1) GeV, respectively. For CKM matrix element entering into the transi- 
tion amplitude, iVf^V^*! = 0.041 is used. The main input parameters of QCD sum rules 
for the form factors are DA's of the S baryon, whose explicit expression are presented 
in the Appendix. Here, we should make the following remark. In the matrix element, 
e"*'^(0|M5^(0)sg(x)(i^(0)|E(p)), besides the functions presented in the Appendix, there appear 
also the functions Af^ , Vf^ and 7^*^ whose explicit forms are unknown for the S baryon. 
Considering the SU (3) / symmetry, we get them from the nucleon DA's. Our calculations 
show that their contribution constitutes only few percent of the final results, so we neglect 
their contribution in the present work. 

Besides these input parameters, there appear also three auxiliary parameters in the sum 
rules, i.e. Borel mass parameter M^, continuum threshold Sq, and the general parameter /3 
arising in the interpolating current of the baryon. These parameters should not effect 
the values of the form factors, so one should obtain working regions of them for which the 
form factors show weak dependence on these parameters. 

Lowering the value of the Borel mass increases the contribution of the higher twist DA's, 
hence requiring that the twist expansion converges leads to a lower limit on the Borel mass; 
on the other hand, increasing the value of the Borel mass increases the contribution of the 
higher states and the continuum. Hence requiring that the contribution of the higher states 
and continuum to the correlation function is less than half the total contribution yields 
an upper bound on the Borel mass. Both of these conditions are met if the Borel mass 
is chosen in the interval 15 GeV^ < < 30 GeV^. The continuum threshold Sq is not 
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totally arbitrary but it is correlated to the energy of the first excited state. Our analysis 
shows that in the region {m^^ + 0.3 GeV)"^ < Sq < {rriY,^ + 0.7 GeVY, the dependences of 
the form factors on this parameter are weak. Finally, to find the working region of (3, the 
dependence of the form factors on cos 6' in the interval —1 < cos6 < 1, where tan^ = /3 is 
considered. Our numerical calculations lead to the working region, —0.5 < cos 6 < 0.7. 
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Figure 1: Analysis of the sum rules for the form factor fi{q'^ 



As an example, in Fig. 1, we depict the dependences of the form factor fi{q^) on 
auxiliary parameters as well as q^. Fig. (la) shows the dependence of this form factor on 
cos^ at the fixed values = 13 GeV^, = 20 GeV^ and sq = (40 ± 1) GeV^. As it is 
seen, there is a stable region in the interval —0.5 < cos^ < 0.7. In Fig (lb), the Borel mass 
dependence of the same form factor is depicted at = 13 GeV^ and the same value of 
the continuum threshold. In the chosen working region of the Borel mass, our predictions 
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change by approximately 5%. From these two figures, it is also seen that our predictions 
are almost independent of the continuum threshold. Finally, in Fig (Ic), we show the 
dependence of the form factor fi on q"^ at two fixed values of /3, and = 20 GeV^ and 
So = 40 GeV^. 

The sum rules predictions are only reliable for the region ^ ml , where in the decay of 
Sfo, the allowed range of extends until [mj]^ — ms)^. To extend the sum rules predictions 
to the whole physical region, the sum rules predictions are fitted to the following function: 

^ m'i ., ' ^ rni.. ' 

fit fit 

The central values of the fit parameters a, b, and mfu are presented in Table [1] This 
table also exhibits values of the form factors at = 0. These errors presented in this 
table are due to the variation of the auxiliary parameters, M^, sq, and /3, as well as the 
errors in the input parameters. Note that for all form factors, the fit mass is always in 
the range nifit = (5.1 — 5.4) GeV. In a vector dominance model, although the double pole 
structure would not be expected, these form factors would have poles at the masses of the 
(axial)vector meson that couples to the transition currents. The observed (axial)vector Bg 
mesons have masses in the range (5.4 — 5.8) GeV. Although the fit mass tends to be slightly 
smaller than the observed masses, considering the uncertainties inherent in the sum rules 
calculations, the results are reasonable. To improve the results one should consider the as 
corrections to the distributions amplitudes and more accurately determine the DA's of S 
baryon. 

Finally we calculate the differential and total decay rate of the T^f, — )■ transition. 
The general form of the differential rate for the rare baryonic weak decay is given by [S3] : 



dV 
ds 



G^a" 



81927r 



m + 3 A(^) 



(20) 



where s = g^/m|^, G = 1.17 x 10 ^ GeV ^ is the Fermi coupling constant and A = A(l, r, s) 
with A(a, b, c) = + 6^ + — 2ab — 2ac — 2bc is the usual triangle function. Here, 

V = ■sjl — is the lepton velocity. The functions 9(s) and A(s) are given as: 



e(s) = 32m^m|^s(l + r - s) 



\D,f + \E,\') 
+ 64m^m|^(l - r - s) Re[DlE3 + D3EI] 
+ QAml y^{Qm^^ - m| s)Re[DlEi] 



+ 64m2m|^v^(^2ms,sRe[D*E3] + (1 - r + s)Re[D*D3 + E^E^] 
+ 32m|^(2m^2 ^ m|^s)|(i _ r + s)m^,y/^ Re[AlA2 + BIB^] 
- ms,(l - r - s) Re[AlB2 + A^Bi] - 2v^fRe[A*5i] + m| sRefA^^a] 



+ 8m|^|4m^(l + r — s) + m 
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+ 8mi <^ 4m' 



A + (1 + r - s)s 



AmJl + r — s) 



(1-r? 



- s']] {\D,f + \E,f) 
+ 8m|^st;=^| - 8m^^sy/^Re[D;E2] + 4(1 - r + s)^/¥Re[DlD2 + EIE2] 
- 4(1 - r - s) Re[DlE2 + D^E^] + m^, 



\2 2 

r — s 



(1^21^+1^21')}, (21) 



A {s) = -8m|yA {\Aif + \Bif + \Dif + \Eif) 



■^niy sv'^X 



\A, 



\Eo 



where r = m|/m| and 







-2m,C,) + {h-gi)Cl 


A2 


= (1 ^ 2) , 




^3 


= Ai(l ^3) , 






= A {gi -gi, 




B2 


= 5i (1 ^ 2) , 






= 5i (1 -> 3) , 




Di 


= (/i - ^1) Cio , 




D2 


= A(l ^2) , 




D3 


= A(1^3) , 




El 


= Di {gi -gi) 




E2 


= E,{1^2) , 




E3 


= El (1 ^ 3) . 





(22) 



(23) 



(24) 

Integrating the differential decay rate over s in the interval, Aml/rri^ < s < (1 — ^JrY, 
we get the total decay rates presented in the Table [2l Finally, to obtain the branching ratios, 
one needs the lifetime of the S;, baryon. Although there is no exact information about the 
lifetime of this baryon, it may be informative to take this lifetime approximately at the same 



order of the 6-baryon admixture, (A;,, H^, S^, f2;,) which is r = 1.391 



+0.039 
-0.038 



X 



10-12 s [52|. The 



results of the branching ratios for different leptons are also presented in Table O It is 
seen that the branching ratio for decays into the electrons or muons are more or less the 
same, while the branching ratio for decay into final states containing r lepton is reduced by 
approximately a factor of four. The order of branching fractions show that these channels 
can be detected at LHC. Comparing the presented results in this work with the results of 
any measurements, one can obtain useful information about the nature of Sf, baryon as well 
as new physics effects beyond the SM. 
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a 


b 


rrifit 




fl 


-0.035 


0.13 


5.1 


0.095 ± 0.017 


/2 


0.026 


-0.081 


5.2 


-0.055 ±0.012 


fs 


0.013 


-0.065 


5.3 


-0.052 ±0.016 


9i 


-0.031 


0.15 


5.3 


0.12 ±0.03 


92 


0.015 


-0.040 


5.3 


-0.025 ± 0.008 


93 


0.012 


-0.047 


5.4 


-0.035 ±0.009 


fl 


1.0 


-1.0 


5.4 


0.0 ±0.0 


fT 


-0.29 


0.42 


5.4 


0.13 ±0.04 


fl 


-0.24 


0.41 


5.4 


0.17 ±0.05 


91 


0.45 


-0.46 


5.4 


-0.010 ±0.003 


T 
92 


0.031 


0.055 


5.4 


0.086 ± 0.024 


T 
93 


-0.011 


-0.18 


5.4 


-0.19 ±0.06 



Table 1: Parameters appearing in the fit function of the form factors, /i, /2, /a, gi, g2, Qs, 
fl: fli fli 9i: 92 9I f^ll theory for — >■ T,i'^£~ and the values of the form factors 
at — 0. In this Table only central values of the parameters are presented. 
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Appendix A 

In this Appendix, the general decomposition of the matrix element, e'*'"^(0|M^(0)(ig(x)s^(0)|S(p)) 
as well as the DA's of E |19] are presented. Considering Lorentz and parity invariances, 
the matrix element can be decomposed into various Lorentz structures as: 

4(0|e'^^X(aia;)4(«2a;)rf^(a3x)|S(p)) 

2 2 

+ V2mE(j^C)„/3(^75S)^ + V3ms(7/.C')„;3(7'^75S)^ + Viml{:i:C)^p{-i^T)^ 

an 



2^2 



+ 



X m 



+ A"^|(;j^75C)Q/3S^ + ^5m|(7^75C)a/3(ia^''x,,S)^ + A"^|(^75C)a/3(;i^S)^ 



2™2 



+ (ri + 



+ r7m|(a^,C),^(a^";i^75S)^ + r8m|(xV^,C),;3(a^''x,75S), 



(Al) 



The calligraphic functions in the above expression do not have definite twists but they 
can be written in terms of the S distribution amplitudes (DA's) with definite and increasing 
twists via the scalar product px and the parameters Oj, i = 1,2, 3. The relationship between 
the calligraphic functions appearing in the above equation and scalar, pseudo-scalar, vector, 
axial vector and tensor DA's for E are given in Tables [3l HJ [5l [6] and [TJ respectively. 



2pxS2 = 5*1 - 5*2 



Table 3: Relations between the calligraphic functions and S scalar DA's. 



Pi = Pi 



2pxV2 = Pl-P2 



Table 4: Relations between the calligraphic functions and S pseudo-scalar DA's. 
Every distribution amplitude F{aipx)= Si, Pi, Vi, Ai, Ti can be represented as: 

F{aipx) = j dxidx2dx3S{xi + X2 + X3 - l)e"*^''^'^'"'F(xj) . {A.2) 

where, with 2 = 1,2 and 3 are longitudinal momentum fractions carried by the partici- 
pating quarks. 
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Vi = Vi 

2pxV2 = V^-V^-V., 

2V3 = 

ApxVi = -2Vi + + ^4 + 2^5 

4pxV5 = Va-V2 

A{pxfV^ = -Vi + V2 + Vz + VA + V^-V(i 



Table 5: Relations between the calligraphic functions and S vector DA's. 

Ai = Ai 

2pxA2 = -Ai + A2-AS 

2A3 = A3 

ApxAi = -2Ai - A3- Ai + 2A5 
ApxA^ = A3- A4 
4{pxfA6 = - + A3 + - +~A^ 



Table 6: Relations between the calligraphic functions and S axial vector DA'; 



The explicit expressions for the S DA's up to twists 6 are given as follows 
Twist-3 distribution amplitudes: 



Vi{Xi) = 120X1X2X30°, Ai(Xi)=0, 
Ti(Xj) = 12OXiX2X303° 

Twist-4 distribution amplitudes: 



Si{xi) 


= 6(x2-xi)x3(e4 + e4°), 


Pl{Xi] 




= 24XiX204 , 


A2[Xi 


V3iXi) 


= 12X3(1 -X3)^t 


A3{Xi 


T2iXi) 


= 24XiX204° , 


T'iixi 


Trixi) 


= 6x3(1 -X3)(e;°-e4°) • 





6(X2 - Xi)x3(^4 - t 



0, 





4 5 



-12X3(X1 - X2)^, 

6x3(1 -X3)(e4° + e4°) 



Twist-5 distribution amplitudes: 

^2(x.) = ^(xi-x2)(e5° + e;°), P2{x.;) = \{x,-X2){il-(^) 

Vj^Xi) = 3(1-X3)V'5, A4(Xi) = 3(xi - X2)V'5 

V^iXi) = 6X30^, A5{Xi) = 0, 

3 
2 

3 



T4{X^) = --i^l+X2)i^','+e,), r5(x,)=6X30;°, 



Tsix.) = 2(^i + ^2)(4'-C 
Finally, twist-6 distribution amplitudes: 



Ve{x^) = 20^, Ae{xi) 



TJxA = 2d)'^. 
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Tl=Ty 




2pxT2 = Ti + T2 


-2T3 






2px% = Ti - T2 


- 2T7 




2px% = -Ti + T 


+ 2Ts 


4:{pxfr6 = 


-- 2T2 - 2T3 - 2T4 - 


f 2T5 + 277 + 2T8 




ApxTj = T7 — 


Ts 


4:{pxfrs 


= -Ti + T2 + T5 - 


- Tg + 2T7 + 2Ts 



Table 7: Relations between the caUigraphic functions and E tensor DA's. 



where, 



0^ = 0^ = ^fe + Ai), 03° = 06° = -e5° = ^(4A3-A2) 

= e4° = ^(8A3-3A2), = -e^o = 1a2 , 

= ^(12A3-5A2), 



/s = (9.4±0.4) X 10"^ GeV^ Ai = -(2.5 ± 0.1) x 10"^ GcV^ 

A2 = (4.4 ± 0.1) X 10-2 GeV^ A3 = (2.0 ± 0.1) x 10"^ GeVl 
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